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ABSTRACT 

We report the discovery of PSR J 1838-0537, a gamma-ray pulsar found through a blind search of data from 
the Fermi Large Area Telescope (LAT). The pulsar has a spin frequency of 6.9 Hz and a frequency derivative 
of -2.2 x 10 -11 Hz s , implying a young characteristic age of 4970 years and a large spin-down power of 
5.9 x 10 36 erg s" 1 . Follow-up observations with radio telescopes detected no pulsations, thus PSR J1838-0537 
appears radio-quiet as viewed from Earth. In September 2009 the pulsar suffered the largest glitch so far seen in 
any gamma-ray-only pulsar, causing a relative increase in spin frequency of about 5.5 x 10~ 6 . After the glitch, 
during a putative recovery period, the timing analysis is complicated by the sparsity of the LAT photon data, 
the weakness of the pulsations, and the reduction in average exposure from a coincidental, contemporaneous 
change in the LAT's sky-survey observing pattern. The pulsar's sky position is coincident with the spatially 
extended TeV source HESS J1841— 055 detected by the High Energy Stereoscopic System (H.E.S.S.). The 
inferred energetics suggest that HESS J1841-055 contains a pulsar wind nebula powered by the pulsar. 
Subject headings: gamma rays: stars - pulsars: individual (PSR J 1838-0537) - ISM: individual objects 
(HESS J1841-055) 



1. INTRODUCTION 

In the last three years, the Large Area Telescope (LAT; 
lAtwood et al.l 120091) aboard the Fermi satellite has proven a 
revolutionary detector of gamma-ray pulsars[3 These objects 
are rapidly rotating, highly magnetized neutron stars, whose 
rotation carries the gamma-ray emitting regions past an ob- 
server's line of sight, creating periodic pulsations. 

With the LAT, pulsars have been detected via three differ- 
ent approaches. First, for pulsars known from radio obser- 
vations, gamma-ray pulsations are detected by coherent fold- 
ing of the LAT photon a rrival times based on the provided 
pulsar ephemerides (e.g.. lAbdo et alJlTolOd : iGuillemot et alJ 
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120121) . Second, radio searches of unassociated gamma-ray 
sources , as listed in the F ermi-LAT Second Source Catalog 
(2FGL: iNolan et aT]l2012l) . have revealed new radio pulsars. 
The ephemerides obtained from radio timing observations 
are in turn used to coherently phase-fo ld the LAT data and 
probe for gamma-ra y pulsations (e.g.: iRansom et alj 1 201 1 ; 
Cognard et all l20llt iCamiloet all 1201 It iKerretalT 12012 : 



RavetalJl2012l) 



Contrasting the previous two strategies, "blind searches" 
are not guided by prior knowledge of the pulsar param- 
eters (see e.g. iChandler et al.l 120011) . Within a year after 
launch of Fermi, pr evious blind searches had impressively 



unvei led 24 pulsars dAbdo et al.l l2009al : ISaz Parkinson et a! 



201C ) using a time-diff erencing technique (lAtwood et al 



20061: iZiegler et al.H2008|). Another two p ulsars were found 



after two years dSaz Parkinson et alj|2012l) . however the de 
tection rate had dropped since then, predominantly because 
of the increasing computational challenge with the longer data 
time span. 

The sensitivity of blind searches for gamma-ray pulsars 
is limited by computing cost, because a grid in the search 
parameter space (for isolated systems typically sky loca- 
tion, frequency /, and spin-down rate /) must be explicitly 
searched. The number of grid points required to discretely 
cover the relevant parameter space i ncreases as a high power 
of the coherent integration time (e.g jBradv et aljfl998l) . This 
problem is analogous to blind searches for continuous grav- 
itational waves (CWs) emitted from r apidly spinning neu- 
tron stars, where hierarchical strategies dSchutz & Papall20Qr ' 
Brady & Creightonl l2000l iKrishnan et al.l 12004 ICutler et al 



2001 " and semi-coherent methods are used to efficiently scan 



wide parameter-space ranges in years of data at finite compu- 
tational resources. 

Here, we present the discovery and key parameters of a 
young, energetic gamma-ray pulsar, PSR J 1838-0537, de- 
tected during a new blind survey e xploiting a novel data- 
analysis method dPletsch et alj|2012h originally developed to 
blindly search for weak CW signals dPletsch & Allen! 
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iPletschl 120 101 1201 ll) . The pulsar is spa tially coincident 
with t he extended source HESS J1841— 055 ( lAharonian et al.l 
120081) . detected at very high energies (VHE; E > 0. 1 TeV). We 
show that the pulsar's energetics make this association plau- 
sible. 

2. DISCOVERY OF PSR J1838-0537 

In recent work dPletsch et al.ll2012h we reported on the dis- 
covery of nine new gamma-ray pulsars in this novel, ongoing 
blind survey using 975 days of LAT data. Unlike those nine 
pulsars, PSR J 1838-0537, being significantly younger, was 
initially detected only during the first year of data (before its 
large glitch). After the glitch, recovering the (weaker) pul- 
sations proved to be a substantial challenge requiring a dedi- 
cated follow-up study presented below (Section[3]l. 

In the new blind-search effort, detecting PSR J 1838-0537, 
we have utilized a hierarchical (multistage) approach. The 
first stage is semi-coherent, where coherent Fourier powers 
computed over a 6-day window are incoherently combined 
by sliding the window over the entire 975 days of data. A 
fundamental new element of the sea rch method is the ex- 
ploitation of a parameter-space metric dPletsch & Allenll2009t 
IPletschl 120101) to build an efficient search grid (in /, / and 
sky position). In a second stage, significant semi-coherent 
candidates are automatically followed up via a fully coherent 
analysis. A third stage further refines coherent pulsar candi- 
dates by using highe r signal harmonics adopting the //-test 
dde Jager et al.lll989l). Further advances incorporated are de- 
scribed in lPletsch et al.l ( 120121) . 

The survey targets unassociated sources from the 2FGL 
catalog with "pulsar-like" properties. Such sources feature 
significantly curved emission spectra and low flux variabil- 
ity over time. Further details on the survey regarding source 
selectio n, data preparation and search parameter-space are 
found in lPletsch et all (120 121) . 

We have detected PSR J 1838-0537 as part of the survey in 
the blind search of 2FGL J 1839.0-0539. This source also has 
gamma-ray counterparts 1F GL J1839.1-0543c i n the Fermi- 
LAT First Source Catalog (lAbdo et al.l l2010al) and OFGL 
J1839 0-0549 in the Fermi Bright Source List ( lAbdo et al.l 
I2009bl) . 

3. FOLLOW-UP ANALYSIS 
3.1. Data Preparation 

For the follow-up study, we prepared a dedicated LAT data 
set extending the blind-search input data to 1 168 days, up to 
2011 October 17. Using the Fermi Science Toolo, we se- 
lected LAT photons whose reconstructed directions lie within 
15° of the pulsar, have energies above 100 MeV, and zenith 
angles < 100°. We included only "Source"-class photons ac- 
cording to the P7_V6 Instrument Response Functions and ex- 
cluded times when the satellite's rocking angle exceeded 52°. 

To improve the signal-to-noise ratio, each photon is as- 
signed a weight meas uring the pr obability that it has origi- 
nated from the pulsar dKerrll201lh . These weights are com- 
puted through a binned likelihood analysis (as done for the 
original blind-search input data) using a spectral model for the 
region including all 2FGL-catalog sources found within 15° 
of PSR J 1838-0537. The pulsar spectrum is modeled as an 
exponentially cutoff power law, dN/dE oc £" r exp (-E /E c ), 
where T is the photon index and E c is the cutoff energy. 

16 http://fermi.gsfc.nasa.gov/ssc/data/analysis/scitools/overview.html 



TABLE 1 
Parameters of PSR J1838-0537 



Parameter Value 



Right ascension, a (J2000.0) 18 h 38 m 56?02(3) 

Declination, 5 (J2000.0) -05°37'09"(2) 

Galactic longitude, I (°) 26.5 

Galactic latitude, b (°) 0.2 

Spin frequency, / (Hz) 6.863015715(4)" 

6.86305339(l) b 

Frequency 1st derivative, / (10~ n Hz s' 1 ) -2.18964(6)" 

-2.2222(l) b 

Frequency 2nd derivative, / (10~ 22 Hz ,v~ 2 ) 5.0(4)" 

17.0(9) b 

Epoch (MJD) 55100.0 

Data span (MJD) 54702 - 55836 

Characteristic age, r c (yr) 4970 

Spin-down power, E (ergs -1 ) 5.9 X 10 36 

Surface magnetic field strength, B s (G) 8.3 X 10 12 

Light-cylinder magnetic field strength, Blc (kG) .... 24.7 

Photon index, T 1.8 ±0.1 

Cutoff energy, E c (GeV) 5.6± 1.4 

Photon flux above 100 MeV (10~ 8 photons cnr 2 s~') 23.9±5.1 

Energy flux above 100 MeV (10 -11 erg cm" 2 s~') .. 17.5 ± 1.8 



NOTE. — Numbers in parentheses are statistical la errors in the last digits. 
a Pre-glitch solution (before MJD 55100). 
b Post-glitch solution (after MJD 55450). 

The extragalactic diffuse emission and the residual back- 
ground are modeled jointly using the iso _p7v6source tem- 
plate, and the Galactic diffuse emission is modeled using the 
gal_2yearp7v6_v0 map cube. For computational feasibility, 
only the 165 000 photons with the highest weights are consid- 
ered (implying a probability-weight threshold of 0.018). 

3.2. Pulsar Timing and Glitch Analysis 

Using the 1168-day data set, the initial pulsar parameters 
are further refined thro ugh a timing-anal ysis procedure. With 
techniques described in lRav etail d201 ll) . pulse Times Of Ar- 
rival (TOAs) are precisely measured and parameters of a tim- 
ing model are fit to these measurements. 

TOAs are obtained from subdividing the follow-up data 
set into 40 segments of about equal length. The initial pul- 
sar parameters are used to fold the photon arrival times and 
produce a set of pulse profiles. By cross-correlating each 
pulse profile with a multi-Gaussian template derived from fit- 
ting the entire data set, the TOAs are det ermined u s ing th e 
unbinned-maxi mum-likelihood m ethod of iRay et al.l (1201 ll) . 
Then Tempo2 dHobbs et al.ll2006l) is used to fit the TOAs to a 
timing model including sky position, frequency and frequency 
derivatives. 

As is typical for a youthful pulsar, PSR J 1838-0537 shows 
irregularities in its spin-frequency evolution that are of two 
types: (i) timing noise, observed as slowly varying, non- 
deterministic fluctuations in frequency, which can be mod- 
eled by including a second frequency derivative in fitting the 
TOAs; (ii) glitch activity, an abrupt spin-frequency increase, 
typically followed by a recovery phase towards the pre-glitch 
rotational state. Thus, the glitch of PSR J 1838-0537 leads 
to rapid loss of phase coherence in the timing analysis if not 
additionally accounted for. 

The spin-parameter changes after the glitch are estimated in 
a dedicated analysis. Fixing the sky position at the pre-glitch 
(before MJD 55100) timing solution, ranges in frequency /, 
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A/(10- 5 Hz) A/ (10- 12 Hz s- 1 ) A/ (l()- 21 Hz s- 2 ) 

FIG. 1, — Pulsar glitch analysis. Color-coded is the weighted //-test statistic for photons lying within a 90-day window that is slid over the entire 
data set with 90% overlap. Fixing the sky position at the pre-glitch (before MJD 55100) solution, for each window scans in //-test over {/,/,/} are done. 
Vertical axes show the time midpoint of each window. Horizontal axes show the offsets from the inferred pre-glitch parameters in / (left), / (middle), and / (right). 



frequency derivative /, and second frequency derivative / are 
scanned on a dense grid ar ound their p re-glitch values, com- 
puting the weighted //-test (Kerr 2011 ) at each grid point us- 
ing photons within a fixed time window. As a balance be- 
tween signal-to-noise ratio and time resolution, we use a win- 
dow size of 90 days, which is just long enough to still produce 
a detectable signal-to-noise ratio. This 90-day window is slid 
over the entire data set with 90% overlap between subsequent 
steps. The results (Figure [TJ indicate that the pulsar experi- 
enced a glitch near MJD 55100. However, immediately after 
the nominal glitch epoch no convincing signal is detected for 
a period of about 300 days. This might be due to the pulsar's 
frequency changing faster during early glitch recovery than 
the {/,/,/} phase model used. In this respect, shorter win- 
dow sizes would possibly help, but these are prohibited by the 
sparse sampling of LAT photons. 

From Figure Q] estimated values for the glitch parameters 
are added to the timing model, serving as input to iterate the 
above timing procedure. Table Q] shows the obtained phase- 
coherent timing solution. The measured post-glitch spin pa- 
rameters imply a relative increase in / of A/// ~ 5.5 x 10 
and a relative decrease in / of -A/// ~ 1.5%. In addition, 
the timing solution requires a relative increase in / of A/// ~ 
2.4, which likely accounts mainly for timing noise. In terms 
of A///, this gl itch is in the top 5% of all pulsar glitches 
recorded to date ( lEspinoza et al.H2.Ql 11), and is the largest ob - 
served for any gamma-ray-only pulsar dDormodv et al.l20"T2h . 
Measuring such strong pulsar glitches is important, as they 
may allow probing the physics of neutron-star interiors (e.g., 
iShapiro & Teukolsky|[T983l; lHaskell et"al]|20T2l) . 

The measured parameters listed in Table Q] characterize the 
pulsar as young and energetic. The implied characteristic age 
and surface magnetic field strength are r = —f/2f = 4970 yr 
andfl s = 3.2x 10 19 (-///V 1 ) 1/2 G= 8.3 x 10 12 G. Assuming 
a neutron-star moment of inertia of / = 10 45 g cm 2 , the pulsar's 
spin-down power is derived as E = -4ir 2 Iff = 5.9 x 10 36 erg 
s" 1 . 

3.3. Pulse Profile and Spectral Parameters 

Figure [2] shows the phase-time diagram and pulse profile 
using the timing solution of Table Q] In the integrated pulse 



profile (weighted pulse phase histogram) statistical errors are 
obtained as C^jwj) 1 ^ 2 , where j runs over all photons in the 
same phase bin and Wj denotes the y'th photon's probability 
weight. Comparing the gamma-ray emission before and af- 
ter the glitch, we found no observable changes in flux, pulse 
profile or spect rum of the pulsar, con sistent with other LAT- 
pulsar glitches (iDormody et al.ll2012l) . 

The integrated pulse profile of PSR J1838-0537 is fitted 
by two Lorentzian lines. The Full-Widths at Half Max- 
ima (FWHM) for the first (PI) and second (P2) peak are 
0.18 ±0.09 and 0.13 ±0.05, respectively. The separation be- 
tween the two is A = 0.24 ± 0.04. Figure [2] also exhibits a 
decrease in the ratio of peaks PI and P2 with in creasing en- 
ergy. Also observed in other LAT pulsars (see e.g. lAbdo et al.l 
1201 Oct iNoutsos et al.l20TTtlAbdo et al.l2010dl) . this is thought 
to be caused by varying gamma-ray emission altitudes and 
curvature radii of th e magnetic field lines as the pulsar rotates 
(lAbdo et al.ll2"0l0dh . 

Observing pulsed emission out to £ max ~ lOGeV im- 
plies transparency to magnetic pair creation. This bounds 
the altitude r of emission (e.g., lAbdo et al.l IMTObh . giving 
r > (E max B n /2.7 GeV) 2/5 p- l/5 R* in flat spacetime (Story 
& Baring, in preparation) for surface fields 10 12 Bi2G. For 
PSR J 1838-0537 one obtains r > 6.1/?*, precluding emission 
near the stellar surface. 

Table[T]lists the best-fit values for T, E c , and the photon and 
energy fluxes of PSR J 1838-0537 as derived from a spectral 
analysis of the region, restricted to photons with pulse phases 
between and 0.5 maximizing the pulsar's contribution. Fur- 
ther details on the spectral analysis of PSR J 1838-0537 and a 
search for unpulsed gamma-ray emission from a putative pul- 
sar wind nebula (P WN) in the off-pulse ph ase interval will be 
reported elsewhere (lAharonian et al. 20121). 

3.4. Pulsation Significance 

Over the 1168 days, the obtained timing solution yields a 
weighted //-test of 486.4, being extremely significant with a 
single-trial false alarm probability of ~ 10~ 92 (Kerr 201 1). 

Figure [3] investigates how this //-test value accumulates 
over time. As expected, the signal increases linearly during 
times where the timing solution describes the pulsar's rota- 
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Fig. 2. — Left: Phase-time diagram, showing the pulse phase for each photon arrival time with the photon probability weight represented by gray-scale 
intensity. Right: The upper plot shows the integrated pulse profile (40 bins per rotation); error bars represent 1<j statistical uncertainties. The five plots below 
show integrated pulse profiles in different energy ranges. For clarity, horizontal axes always show two rotations. 



tional behavior well (before the glitch at MJD 55100 and af- 
ter MJD 55450). However, comparing the two different slopes 
from linear fits made during both time intervals, the slope dur- 
ing the post-glitch interval is smaller by about 27%. 

This effect is traced back to a reduction in exposure result- 
ing from a change in /•errai-LAT's sky-survey observing pat- 
tern on 2009 September 3 (MJD 55077). This change has 
decreased the average exposure at the pulsar's sky position 
by about 24% (FigureO middle panel). Hence, the exposure 
reduction explains (within 3%) the smaller rate in //-test in- 
crease over the later data (the //-test is linear in the number 
of photons). This is further confirmed by observing no signif- 
icant variation in weighted photon flux over time (Figure [3] 
lower panel). 

4. MULTIWAVELENGTH OBSERVATIONS 

4.1. Radio and X-ray Counterpart Searches 

To determine whether PSR J 1838-0537 is also visible as a 
radio pulsar, a 1 .7-hr observation has been conducted on 201 1 
May 27 with the Green Bank Telescope (GBT) at 2 GHz. 
We also analyzed a 1-hr archival GBT observ ation from on 
2009 July 5 at 0.8 GHz (iRansom et al.l 1201 ll) . Configura- 
ti on details of both ra dio observations are found in Table 5 
of iPletsch et al.l (120121) . Folding the data using the timing so- 
lution of Table Q] and only searching in dispersion measure 
revealed no radio pulsations. 



As in IPletsch et al.l (120121) . upper limits on the radio flux 
density are derived with the modified radiometer equation 
(lLorimer & Kramerl 120051) . The sensitivity loss due to the 
telescope pointing-direction offset from the pulsar position 
(0.5' and 5.7', respectively) is app ropriately accounted for, 
analogously to IPletsch et al.l (120121) . leading to radio flux up- 
per limits of 9 /xJy for the 2 GHz observation, and 82 /iJy 
for the 0.8 GHz observation. These upper limits are com- 
parab le with those ob t ained for other LAT -discovered pul- 
sars dRav et all 120111: IPletsch et all 120121) suggesting that 
PSR J1838-0537 is radio-quiet. 

In searching for X-ray counterparts, two archival Swift 
observations (4.5 ks total exposure) cover the position of 
PSR J 1838-0537, but no counterpart is observed. Using a 
41.1-ks Suzaku observation, we find a hint of X-ray emission 
at the pulsar location. However, the Suzaku large positional 
uncertainty (<~ 30") and the low confidence (~ 3er) of the 
counterpart candidate preclude a firm X-ray detection of the 
pulsar. The unabsorbed flu x upper limit i n the .3-10 keV en- 
ergy range computed as in iMarelli et all d20TTh is 5.0 x 10" 14 
erg cm~ 2 s -1 , simila r to other LAT pulsars (IMarelli et alJl201 it 
IPletsch et alJLMl) . 

4.2. A Likely Pulsar Wind Nebula Association 

The sky position of PSR J 1838-0537 coincides with the 
unidentified VHE source HESS J1841-055 discovered dur- 
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Time (10 4 MJD) 

FIG . 3 . — Upper panel: Weighted //-test statistic as a function of time using 
the timing solution of Table [l] (dotted curve). Black solid lines represent 
separate linear fits before the nominal glitch (before MJD 55100), and after 
recovery of pulsations (from MJD 55450) where the measured slope is 
smaller by 27%. Middle panel: Exposure for the pulsar's sky location in time 
bins of size 6.7 days. The average exposure (horizontal dashed line) after the 
LAT's observing-pattern change at MJD 55077 (vertical dashed-dotted line) 
is smaller by 24%. The exposure oscillation is at the 53.4-day precession 
period of the spacecraft orbit. Lower panel: Sum of photon probability 
weights per time bin divided by exposure. Binning in time is identical to the 
middle panel and error bars represent \a statistical uncertainties. 



ing t he H.E.S.S. Galactic Plane Survey dAharonian et alj 
120081) . Although its extended TeV morphology suggests that 
HESS J1841-055 is composed of multiple sources, limited 
exposure has hitherto prohibited statistically significant con- 
firmatio n. None of the previou sly known counterparts pro- 
posed in lAharonian et alj (12008b is capable of fully account- 
ing for the VHE emission alone. 

PSR J1838-0537 is sufficiently energetic to power a PWN, 
further supporting this association. With distance d to the 
pulsar, the spin-down flux at Earth is E/d 2 . If assuming 
a pseudo-distanc^] based on the observed correlation be- 
tween t he gamma-ray luminosity (0.1- lOOGeV) and E (fol- 
lowing ISaz Parkinson et al.ll2.010t) of 2kpc this yields a large 
spin-down flux E/d 2 ~ \ .5/d 2 _ kpc x 10 36 erg s" 1 kpc" 2 (where 
dikvc = fif/2kpc), wh ich is comparable to similar systems 
dCarrigan et alj200g) . 

The integral energy flux of HESS J1841— 055 at the detec- 
tor over t he range 0.5 -80 T eV is G V he ~ 5.8 x 10" 11 erg 
cm" 2 s" 1 dSguera et al.l 120091) . Assuming isotropic emission 
this implies a conversion efficiency rj = L 1 jE = Q.5d 2 kpc % 
(0.5 -80 TeV), similar to ot her suggested pulsar/PWN asso- 
ciations dHessel s et al. 2 0081) . 

4.3. A Nearby Candidate Supernova Remnant 

Because the pulsar is very young, we also consider the 
possible association with the nearby candidate supernova 

17 This pseudo-distance is subject to various caveats that trans- 
late to considerable unce rtainties in this estimate of factors of a few 
( Saz Parkinson et al. 2010). 



remnant (SNR) G26.6-0.1. This diffuse X-ra y source was 
detec ted in the ASCA Galactic Plane Survey dBamba et alj 
12003b . T he best-fit ab s orptio n column density of the X-ray 
spectrum iBamba et alj d2003i) yields a distance estimate of 
1.3 kpc. 

The SNR candidate has an angular size of about 12' 
(FWHM), and its approximate geometric center lies 20' away 
from the pulsar's sky position. We note that the probabili ty of 
chance superposition is also not negligible (Kaspi 199M). To 
assess the SNR association, we assume the 1.3 kpc distance 
and an age of 5 kyr, and consider the corresponding pulsar 
transverse velocity obtained as ~ 1500 (d/ 1.3 kpc ) km s -1 . 
This is about 50% faster tha n the highest velocity directly 
measured for a neutron star dChatteriee et alj 120051) . How- 
ever, the distance is quite uncertain and could be significantly 
smaller, which would reduce the resulting transverse veloc- 
ity. Thus the SNR association appears unlikely, but cannot be 
fully excluded. 

5. CONCLUSIONS 

We report the discovery and follow-up study of 
PSR J 1838-0537 found in a new blind-search effort us- 
ing Fermi-LKT data. The inferred parameters distinguish 
the pulsar as young and energetic. The characteristic age, 
T = 4970yr, makes PSR J1838-0537 the second-youngest 
pulsar ever found in a blind search of gamma-r ay data, after 
PSR J 1023-5746 dSaz Parkinson et alj l20lob . for which 
t = 4610yr. Considering all rotation-powered pulsars known, 
PSR J 1838-0537 resides among the top 1% of systems with 
largest E. In deep radio observations no pulsations have been 
detected, suggesting that PSR J 1838-0537 is radio-quiet, 
only detectable through its gamma-ray pulsations. Moreover, 
the pulsar experienced the largest glitch yet observed in any 
gamma-ray-only pulsar. 

The pulsar's sky location coincides with the extended VHE 
source HESS J1841-055. We have shown that the pulsar's 
energetics are likely to power a PWN producing part of the 
H.E.S.S. -detected TeV emission. This clearly motivates a fur- 
ther analysis of the region at TeV energies to confirm such 
a PWN association. With the LAT, an analysis for off-pulse 
emission of PSR J 1838-0537 will help studying a potentia l 
PWN and the ambient medium (e.g. jAckermann et alj|201 lb . 

PSR J 1838-0537, the tenth object found in this new sur- 
vey, brings t he number of ga mma-ray pulsars undetected in 
radio to 33 dRav et al.l 120121) . Together, these represent a 
significant increase of the gamma-ray-only pulsar population 
and will be an impor tant contribution to the second Fermi- 
LAT Pulsar Catalog dAbdo et al.1 120121). In add ition to the 
improved search techniques dPletsch et al.l 120 12b . the com- 
putational resources for the survey have recently been sig- 
nificantly enhanced by using the volunteer computing sys- 
tem Einstein® Hom^^. The unprecedented search sensitivity 
from the combination of these advances warrants optimism 
for further gamma-ray pulsar discoveries. 
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